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We carried out magnetotransport experiments under hydrostatic pressure in order to study the 
nature of the metal- insulator transition in BaVS3. Scaling relations for p(T,H,p) are established 
and the pressure dependence of the spin gap is determined. Our new results, in conjunction with a 
re-analysis of earlier specific heat and susceptibility data, demonstrate that the transition is weakly 
second order. The nature of the phase diagram in the T-p-H space is discussed. 
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There has been great progress in understanding metal- 
to-insulator transitions which are accompanied by mag- 
netic ordering and/or structural distortions [|l). In con- 
trast, "pure" Mott transitions which do not involve any 
apparent change of symmetry, remain little understood. 
The suppression of magnetic order is usually due to frus- 
tration, or possibly orbital fluctuations. Differences in 
lattice structure and the subspace of the relevant orbitals 
make each of these systems unique. 

The recent availability of single crystal specimens of 
BaVS3 has brought fresh insights |§||], but understand- 
ing the character of its Mott transition, and the nature of 
the low-temperature phases, remains a challenging prob- 
lem. Under atmospheric pressure, the material under- 
goes a metal-to-insulator transition at Tmi = 69K from 
a T > Tmi bad metallic state to a T < Tmi insulator. 
There is a further phase transition at Tx = 30K (0,0, 
and it is now understood that long-period magnetic or- 
der sets in here H. However, in this paper we will be 
concerned with the intermediate (Tx < T < Tmi) in- 
sulating phase (we will call this the insulating phase) 
which is isostructural with the metallic phase, and has 
no magnetic long range order |3|,[ij. It is a singlet insu- 
lator gj] characterized by a singlet-triplet gap. Various 
measurements give varying estimates for the spin gap: 
NMR/NQR (J] suggested A s = 22meV, but according 
to inelastic neutron scattering studies As ~ lOmeV is 
more likely |9[|. In any case, it is remarkable that a spin 
gap As > ksT appears in a 3D system, without the 
formation of any static pattern of bonds. 

Our paper has two main points. First, we use mag- 
netoresistivity measurements to determine the pressure 
dependence of the spin gap up to p — 15kbar. Here we 
closely follow Booth et al. [10] who did a similar analysis 
for p — lbar. However, we managed to improve their 
procedure by dropping some of the simplifying assump- 
tions, and derive corresponding statements from the data 
instead. - Second, we use our fresh data in conjunction 



with previous susceptibility g] and specific heat Q| data 
to discuss the properties of the thermodynamic potential 
G(T,p,H), and determine the character of the metal- 
insulator transition. In contrast to previous claims [[12| . 
we conclude that it is a genuine second order phase tran- 
sition. 

Single crystals of BaVS3 were grown by Tellurium flux 
method. The samples used in this work come from the 
same batch as those investigated in Ref. J2|. The re- 
sistivity measurements under pressure were performed 
in a non-magnetic copper-beryllium cell using kerosene 
as pressure medium. The magnetoresistance (MR) was 
studied in two ways: by magnetic field sweeps at various 
temperatures and by measuring the temperature depen- 
dence of the resistivity in zero and in high magnetic field 
(-ffmax = 120kG). By simultaneous application of carbon 
glass and capacitive thermometers the overall uncertainty 
of temperature was reduced below ±0.05K. 

Figure 111 shows representative p(T) curves measured 
at various pressures both in H = and H = 120kG 
magnetic field. In accordance with earlier observations 
P,|l2| , the dominant effect is the pressure induced reduc- 
tion of the transition temperature. The influence of the 
applied magnetic field is comparatively weak. It is an 
excellent approximation that for all pressures the resis- 
tivity depends on the magnetic field only through the 
field-induced shift of Tmi; the resistivity in H max is the 
shifted counterpart of the zero-field one. It seems a plau- 
sible assumption that any smaller field has the same effect 
on the resistivity, i.e., 



p(T,H,p) = p((T -T MI (H,p),p) . 



(1) 



Tmi itself can be determined with a high accuracy by the 
logarithmic derivative of the resistivity, d(lnp)/d(l/T), 
as shown in the lower panel of Fig. pi 

The field dependence of the transition temperature is 
weak, but it is increasing with increasing hydrostatic 



pressure: 
pressure 



ATMi(-ffmax) grows from 0.35K at ambient 
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FIG. 1. Temperature dependence of the normalized resistiv- 
ity (upper panel) and its logarithmic derivative (lower panel) 
at various pressures. The sharp peaks define the transition 
temperatures at p — 1 bar, 7. 5kbar, and 1 5kbar. 

As expected from the data shown in Fig. [lj the magne- 
toresistance becomes large in the vicinity of the metal- 
insulator transition. Field sweeps performed at constant 
temperatures revealed a quadratic magnetoresistance for 
each pressure, 



Ap(T,H,p)^H 2 



(2) 



Figure @ shows representative curves for p = lbar, and 
p = 15kbar. Note that far above the transition the mag- 
netoresistance is strongly reduced (upper panels). In 
fact, for T > 2Tmi, no effect can be observed, confirming 
that the dominant term in the magnetoresistance is due 
to the field dependence of Tmi and any other contribution 
to Ap is below the detection limit. 

Taken in conjunction with Eqn. (|lj), Eqn. (g) shows 
that at the field strengths available to us, ATmi(H) is 
purely quadratic: 

AT M i(p,H) = -a(p)H 2 . (3) 

In the sense of Eqn. (fy), at any given pressure Ap(H) 
can be expanded in terms in ATmi(H), or effectively in 
powers of H . Up to second order in ATmi 



Ap(T,p,H)= - 



dp(T, P ) 



dT 



AT m (p,H) 



H=0 



i d 2 p(T, P ) 

2 



a{p) 



dT 2 
dp(T,p) 



AT M i(p,H) 2 = 



H=0 



dT 



H 2 + -a 2 {p) 



d 2 p(T,p) 



H=0 



dT 2 



H\ 



(4) 



(5) 



Though ATmi(H) contains no fourth-order shift, the oc 
H A term appearing in Eqn. (ra) is not negligible because 
p(T) is a sharply varying function of T near Tmi- 




H (10 4 xGauss) 
FIG. 2. Magnetic field dependence of the resistivity at var- 
ious temperatures at p — lbar, and p — 15kbar. The dashed 
curves indicate the quadratic fit to the data. Note the different 
scales in the upper panel. 

The consistency of our scheme can be verified, since 
the quantities on either side of Eqn. (||) can be measured 
independently. In Fig. || full squares show the directly 
measured Ap(T,H max )/p(T), as obtained from the dif- 
ference of the p(T) curves measured in 120kG, and in 
zero-field, resp. This set of data is to be compared to the 
curves calculated from (|J) using the partial derivatives 
of resistivities measured at zero field, and inserting the 
value of ATMi(-ffmax) determined experimentally (Fig. [I], 
upper panel). There is no free parameter in the calcu- 
lation. The results of expanding to order (ATmi) 2 are 
shown in Fig. |3j by open circles. The agreement is good. 
It is clear that the second term of the expansion is neces- 
sary to describe the results near Tmi, and that it is also 
sufficient. 

The first order expansion (shown by dashed line in 
Fig. ||) gives already a good approximation except in the 
closest vicinity of the transition. In the range where the 
first order approximation is valid, 



Ap(H,p) 



dp(T,p) 



dT 



AT M1 (H,p) = -(3(p)H 2 (6) 



H=0 



i.e., the magnetoresistance curves determined by field 
sweeps scale to the magnetic field dependence of the tran- 
sition temperature. The magnetoresistance has been de- 
termined at several temperatures at each pressure; the 
corresponding values at iJ max are shown in Fig. S by 
open triangles. This reasserts that the field dependence 
of the transition temperature Tmi(H) is quadratic in H 
for all pressures. 

Eqn. (|3|) can be brought to the dimensionless form 



ATmi(H) (gSp B H 

= = — T — = 



T 



MI 



V k B Ti 



MI 



(7) 



H=0 



where 5=1/2 and g = 2. The sensitivity of the transi- 
tion to the applied field is characterized by the constant 



7. The magnetoresistance results yield 7 
pendently of pressure. 



0.45, inde- 
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p=1 bar 
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FIG. 3. Analysis of the magnetoresistance at various pres- 
sures. Full symbols: Ap(T, -ff max )/p(T) measured by thermal 
sweeps in presence of magnetic field, triangles: results of mea- 
surements by magnetic field sweeps at selected temperatures, 
dashed lines: first order calculation according to Eqn. M); 
open circles: second order calculation from the zero field p(T) 
curves according to Eqn. M). 

We note that an analogous relationship describes the 
suppression of the spin-Peierls transition by a magnetic 
field, with very similar values of 7. On the theoretical 
side, Bulaevskii et al. [|13| predicted 7 = 0.44, while Cross 
Jl4| found 7 = 0.38. Later experiments established that 
7 = 0.41±0.05 is universal both for organic p5[ and inor- 
ganic pq ] compounds. We may hope useful insights from 
such parallels, though the spin-Peierls transition involves 
a static distortion which is not the case with BaVS3. 
Since the magnetic field acts by making an antiparallcl 
alignment of magnetic moments less favourable, we guess 
that the nature of short-range correlations in BaVS3 is 
similar as in spin-Peierls systems. 

Another dimcnsionless form of (0) is obtained by in- 
troducing the pressure dependent spin gap Ag(p): 



ATmi(H) 



Tj 



MI 



( gp B SH 
V A s (p) 



(8) 



assuming that it is related to a characteristic magnetic 
field H c via the corresponding Zeeman energy. We de- 
rived the coefficient a = 1.0 ± 0.1 by accepting the am- 
bient pressure value of the spin gap A sip — latm) = 
lOmeV || . H c can be also understood as the critical mag- 
netic field which would suppress the metal-to-insulator 
transition completely but one has to bear in mind that 
only the weak-field behavior is known (with our assump- 
tion H c — 1700kG). Analyzing their p — latm mag- 
netoresistivity data, Booth et al. found H c = 2600kG 
(175K) [|0|. Let us note, though, that Booth et al. pos- 
tulated a specific (elliptical) form of the phase boundary 
in the T—H plane, and thus in their scheme, a is fixed by 
geometry rather than measured. 



Equation (g) allows to determine the pressure depen- 
dence the spin gap A sip) directly from the measured 
Tmi shifts (Fig. Q). They can also be derived from mag- 
netoresistance measurements performed at constant tem- 
peratures using Eqs. (||) and (§h which relate a to Pip). 
The results are shown in Figure where the pressure de- 
pendence of the transition temperature is also presented. 
We found that within the experimental resolution the 

' "1 the transition 
bT M i(p)- 
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10 

p (kbar) 
FIG. 4. The pressure dependence of the spin gap (di- 
amonds, left-hand scale) scales with that of Tmi (circles, 
right-hand scale). The spin gap determined from the shift of 
the transition temperature and from the magnetoresistance is 
shown by open and full diamonds, respectively. 

The transition temperature Tmi is marked as a peak 
in the d{lnp)/d(\/T) vs T plots but it has been long 
disputed whether it belongs to a true thermodynamic 
transition. The early suggestion of a weakly first or- 
der transition |6| was discarded in favour of a smooth 
conductor- to- insulator crossover ||12| . This might have 
been either of the kind of the supercritical behavior ob- 
served near the critical point of the X^-zCr^O^ system 
Q , or like the behavior of Kondo insulators |17j . How- 
ever, we argue against this interpretation, and show that 
the transition is a genuine phase transition. As to the 
first alternative, the nearness to a critical-point situa- 
tion should be an accidental feature and we should ex- 
pect that changing the parameters drives either towards a 
first-order phase transition, or well away from the critical 
point. However, the resistivity curves show that essen- 
tially the same transition is taking place at all pressures 
up to p cr = 20kbar |pL|l8|, and/or in accessible magnetic 
fields. - As to the second possibility, the crossover from 
the low-temperature insulator to the high-temperature 
conductor in Kondo lattices can be generically continu- 
ous, because the insulating ground state does not violate 
Luttinger's theorem. Actually, this possibility cannot be 
trivially refuted for BaVS3, because the unit cell con- 
tains two V sites, i.e., an even number of 3d electrons. 
However, we explicitly show below that there is a phase 
transition at Tmi- 

Unfortunately, at p > lbar, only resistivity data are 



known, thus the following discussion has to be restricted 
to the ambient pressure case. The lowest-order quan- 
tities which clearly show non-analytic behavior are the 
temperature derivatives of the susceptibility components: 
dxc/dT and dxa/dT have large jumps at Tmi |p|,p^[- If 
it were only for these, d X /dT = -{d 3 G/dTdH T ) be- 
ing a third derivative of the thermodynamic potential 
G(T,p, H) 1 one might have suspected that the transition 
is of third order. We should, however, seek to relate 
A(d X /dT) = (dxi/dT) - (dxu/dT) to other thermody- 
namic quantities by general reasoning [E0| . 

The thermodynamic potential G(T,p, H) must be con- 
tinuous across the transition: f(T,p, H) — Gi(T,p, H\ 
Gm{T,p, H) = along the phase boundary given by 
T M1 (H) = T&j - a H 2 , where T^ = T M1 (p = lbar), 
and ciq = a(p — lbar). Expanding / to fourth order 
in H, and considering that the entropy and the linear 
susceptibility are continuous across the transition, and 
that neither of the phases has spontaneous magnetiza- 
tion Mj(T,H = 0) = M M (T,H = 0) = 0, we are left 
with the relationship 



A 



dx 
dT 



no 
T' 



AC 



1 



12a 



Ax 



(3) 



(9) 



The discontinuity of dx/dT has to be balanced by a com- 
bination of the discontinuities of the specific heat C, and 
the non-linear susceptibility x ■ 

We have measured A(dx/dT) [g|, and clq. AC and 
Ay( 3 ' have to be taken from published data. Imai et al. 
|yj find a strong anomaly at ~69K in the electronic spe- 
cific heat which they chose to identify as a sharp peak 
rather than as a discontinuity but this interpretation is 
not compelling. We have made a new plot using their 
published data points and find that they allow that part 
of the peak height is made up by a discontinuity (Fig. a, 
left). Our fit with (C\ — Cm)/Cm = 1.4 is certainly some- 
what arbitrary, but the correct value cannot be very dif- 
ferent. 

Though the so obtained AC is of the order of magni- 
tude required by (||), a fraction of A(dx/dT) has to be 
matched by Ax^ ■ High- field magnetization data have 
been published by Booth et al. [fill. The authors empha- 
size that the magnetization curves are essentially linear, 
but in fact the data shown in their Fig. 2 are compati- 
ble with a size of the discontinuity Ax^ 3 -* < 0, which is 
sufficient to satisfy (g). 

The character of the transition is described by Eqn. 
(||) which contains both second and third derivatives of 
G. For this reason, one might say that the continuous 
phase transition is "weakly second order" . In any case, 
the existence of a surface of continuous phase transitions 
in the T-p-H space strongly suggests that there is a dis- 
tinct non-magnetic insulating phase which differs from 
the metallic phase not in transport properties only, but 
also in the sense of possessing a hidden order. 
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FIG. 5. Comparison of the the specific heat (based on 
taken from, \lgj) and temperature derivative of the c-axis sus- 
ceptibility (based on \Aj). The lines are guide to the eye. 

In conclusion, magnetotransport measurements under 
pressure were performed on single crystals of BaVS3 in 
order to study the nature of the metal-insulator tran- 
sition. We determined the pressure dependence of the 
spin gap in the insulating phase and showed that As 
scales with the transiton temperature. We discussed the 
nature of the phase diagram and pointed out that the 
metal-insulator transition is not a smooth crossover but 
a genuine phase transition. 
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